Detection of high-energy gamma rays from winter thunderclouds 
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Abstract 

A report is made on a comprehensive observation of a burst-like 7-ray emission from thunder- 
clouds on the Sea of Japan, during strong thunderstorms on 2007 January 6. The detected emission, 
lasting for ~40 seconds, preceded cloud-to-ground lightning discharges. The burst spectrum, ex- 
tending to 10 MeV, can be interpreted as consisting of bremsstrahlung photons originating from 
relativistic electrons. This ground-based observation provides first clear evidence that strong elec- 
tric fields in thunderclouds can continuously accelerate electrons beyond 10 MeV prior to lightning 
discharges. 

PACS numbers: 82.33.Xj, 92.60.Pw, 93.85.-q 
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INTRODUCTION 

As first suggested by Wilson [l| , intense electric fields in thunderclouds and lightning are 
considered to accelerate charged particles to relativistic energies, which in turn will produce 
Bremsstrahlung photons. Thus, detailed investigation of such photons is expected to provide 
a valuable key to the particle acceleration process in strong electric fields. 

Thunderclouds are found to have electric fields reaching 100 kVm -1 |2|, or sometimes 
even > 400 kVm -1 [3|. Usually, frequent collisions in the dense atmosphere will not allow 
electrons to be accelerated even under such a strong electric field. However, seed electrons 
more energetic than ~ 100 keV, if once produced by, e.g., cosmic rays, can gain energies 
from the electric field fast enough to overcome the collisional energy losses. Thus, such 
electrons are accelerated and multiplied in the atmosphere, because their mean free paths 
get progressively longer, up to ~ 50 m, as they acquire higher energies. This is the concept 
of runaway electron avalanche mechanism |4|. When these conditions are fulfilled, electrons 
can be accelerated in thunderclouds to several tens MeV energies or higher, and then produce 
Bremsstrahlung 7-rays. 

Actually, sudden non-thermal x-/7-ray emissions from upper atmosphere, called terres- 
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trial 7-ray flashes, have been observed at equatorial latitudes by near Earth satellites 
The spectra of the flashes are roughly expressed by power-law functions with the photon 
indices of ~ —1 [7|], some of them extending up to 10 — 20 MeV. Now they have been 
interpreted as arising from high-altitude electrical discharges above thunderstorms. In addi- 
tion, ground-based observatories have also detected similar photon bursts from natural js, l| 



and rocket-triggered 
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ll| lightning, revealing that the bursts are mostly associated with 
lightning discharge processes. The burst spectra in the rocket triggered lightning mostly 
extend up to ~ 250 keV, or to 10 MeV [12 



on rare occasions. These bursts including the 
terrestrial 7-ray flashes have typical duration of milliseconds or less. 

In addition to those short-duration bursts, more prolonged radiation from thunder ac- 
tivity, lasting for a few tens seconds to minutes, have been observed by detectors at high 
mountains during summer seasons [13|, [14|. Also radiation- monitoring posts, arranged in 



nuclear power plants in the coastal area of Sea of Japan, have frequently detected such 



prolonged events associated with winter thunderstorms 



15]. However, compared with the 



short duration events, the longer duration ones have remained much less understood, due 
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to inadequate information on, e.g., the particle species, their arrival direction and energy 
spectrum. Here we show one observation revealing those essential features of the prolonged 
burst. 

EXPERIMENT 

Kashiwazaki-Kariwa nuclear power plant is located in the coastal area of Sea of Japan 
in Niigata prefecture, where thunder activity is very high in November through January. 
Installed at the rooftop of a building in this power plant, our new automated radiation de- 
tection system has been continuously and successfully operated since 2006 December 22. The 
system consists of two independent and complementary parts, Detector-A and Detector-B, 
placed with a separation of ~ 10 m. Designed after the Hard X-ray Detector experiment [3] 
on board the 5th Japanese cosmic x-ray satellite Suzaku, Detector-A consists of two iden- 
tical main detectors, each having a structure shown in Fig. [H it comprises a cylindrical 
Nal scintillator operated in 40 keV — 3 MeV, surrounded by a well-type BGO (Bi4Ge30i2) 
scintillator shield counter of which the average thickness is around 2.0 cm. Placed above the 
two main counters is a 5 mm thick plastic scintillator. Detector-A records the energy deposit 
of each event, and its arrival time with a time resolution of 10 /xs. As an important feature 
of Detector-A, the surrounding BGO shields prevent low-energy environmental radiations 
from impinging on the Nal scintillator. Furthermore, using the signals of the BGO shields, 
we can roughly know where an incident particle comes from. When the Nal scintillator is 
triggered by an incident particle and the Nal signals does not coincide with the BGO ones, 
the incident particle would generally arrive from the vertical direction neither the horizontal 
one nor the ground since the BGO shields firmly cover side and bottom of each Nal scin- 
tillator. Thus, using the BGO signals in anti-coincidence, we can reduce background of the 
Nal scintillator, and give it broad directional sensitivity towards the vertical direction. 

Aiming at higher energy ranges, Detector-B utilizes Nal and Csl scintillation counters 
observing in 40 keV — 10 MeV and 600 keV — 80 MeV, respectively. With a spherical 
shape, both work as omni-directional monitors without anti-coincidence shields. Each of 
them accumulates pulse-height spectra in 1024 channels every 6 seconds, and record broad- 
band pulse counts every second. 

In addition to Detectors A and B, we measure environmental light and electric field. Three 
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light sensors are used to measure optical intensity in different directions, each consisting of a 
Si photo-diode (HAMAMATSU S1226-8BK) and a handmade analog circuit. Output signals 
of each light sensor are fed to a 12 bit ADC, and recorded as ADC values every 0.1 sec. The 
electric field sensor is a commercial product (BOLTEK EFM-100), and its analog output is 
collected by a 12 bit ADC every 1 sec and converted to the electric field strength between 
± 10 kVm -1 with a resolution of 5 Vm" 1 . Details of our system are reported by Enoto 
and will be published elsewhere. 



RESULTS 



On 2007 January 6, the coastal area of Sea of Japan was covered with one of the strongest 
thunderstorms in this winter. Figure [2] shows count rate histories of the BGO and Nal 
scintillators of Detector-A, and of the Nal and Csl counters of Detector-B, obtained from 
21:10 to 22:10 UT on this day (6:10 to 7:10 local standard time on January 7). Superposed 
on gradual background changes due to fall-outs of radioactive Radon, we have detected a 
remarkable count enhancement at about 21:43, lasting for ~ 1 min. This event, or "burst", 
is thus highly significant in all the inorganic scintillators of our system. 

Figure [3] gives details of the count-rate histories between 21:40 and 21:50 UT. Records 
of light and electric fields are also shown. The signal in Detector-A is most prominent in 
the > 3 MeV range covered by the last overflow bin [Fig . EJ^a)], where the excess counts 
in a 36 sec of 21:43:09 — 21:43:45 UT, calculated by subtracting the background from the 
total counts during the burst, become 129 ± 17 and 43.6 ± 3.5, without and with the BGO 
anticoincidence, respectively. The Nal and Csl counters of Detector-B recorded 3 — 10 MeV 
excess counts of 74.2 ± 9.6 and 53.3 ± 8.4, respectively [Fig. 13(b)]. Therefore, the burst 
intensity from the two Detectors (without anticoincidence for A) agrees with one another, 
since the values from Detector-A sum the two counters. 

At 21:44:18 UT, or 70 seconds after the burst onset, the optical data rapidly increased 
for about 1 second [Fig. Efc)], with a rise time less than the sampling time which is 100 ms. 
In coincidence with this flash, the electric filed quickly changed its polarity from positive 
values (> 5 kVm -1 ) during the burst to negative (< —10 kVm -1 ), although the signal is 
saturated at ±10 kVm -1 due to the limited range of the electric field sensor. Thus, the 
first flash can be interpreted as a lightning discharge; so are the subsequent four flashes. 
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The burst had already ceased by the first lightning, and no increase accompanied any of the 
total five discharges. Thus, the observed event appears to be associated with thunderclouds 
themselves, rather than to the lightning discharges. 

During the burst, the plastic scintillator of Detector-A [Fig. 13(e)], with its lower threshold 
set at 1 MeV, showed no significant excess signals. Quantitatively, the number of events 
detected by the plastic scintillator with an area of 464 cm 2 in the 36 sec is 185 ± 14, which 
is consistent with the number, 191, predicted by its longer-term average count. The lack of 
excess signals in the plastic scintillator, together with the long-burst duration, excludes the 
burst being due to electrical noise. Furthermore, we can conclude that the burst is dominated 
by 7-rays, since the thin plastic scintillator has high detection efficiency for charged particles 
but is nearly transparent to 7-rays. Another important information, derived from the data, 
is that the photons came from a sky, probably covered with thunderclouds, rather than 
from the ground. Actually, excess count ratio between Nal [Figj2(b)] and BGO [Figf5](a)] of 
Detector-A, ~ 0.18 ± 0.03, is significantly higher than their normal count ratio, 0.08, which 
is due mostly to environmental radioactivity. 

Figure H] shows pulse-height spectra of the burst from Detector-A and B, accumulated 
for the same 36 sec; we subtracted background averaged over 21:30 — 21:40 UT and 21:45 
— 21:55 UT. Thus, the excess counts in both systems exhibit very hard continuum spectra. 
The spectra of Detector-B, in particular, extend up to 10 MeV. Supposing that the 7-rays 
came from the thunderclouds, these spectra imply that electrons were accelerated in the 
thunderclouds at least to 10 MeV, or to higher energies in order to efficiently produce the 
10 MeV photons. 

Based on results from an experiment and from a Monte Carlo simulation [isj], we may 
assume that the photon number spectrum arriving at the top of our system has a power-law 
form as a x E@, where a is normalization in units of photons MeV _1 cm -2 , /3 is the photon 
index, and E is the photon energy in MeV. In this work, the data points in 600 keV — 10 
MeV where both Nal and Csl counters can cover are selected [Fig. IU(b)]. Then, convolving 
this model with the detector responses, and fitting it to the observed spectra (but excluding 
data points with too large errors), we obtained a = 2.63 ± 0.92 and (3 = —1.58 ± 0.23 
(x 2 /d.o.f. = 3.3/3) from Nal of Detector-B, while a = 1.93 ± 0.42 and p = -1.69 ± 0.15 
(x 2 /d.o.f. = 5.3/6) from Csl. Thus, the two sets of results are consistent. Weighted means of 
the two estimates become a = 2.04±0.38 and (3 = — 1.66±0.13. The spectra from Detector-A 
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are generally consistent with this. In reality, the observed photons must be produced at some 
distance, and propagate in the atmosphere while they are attenuated mainly via Compton 
scattering. When corrected for the mild energy dependence of the Compton cross section, 
the photon index softens to —2.01 ± 0.14, —2.44 ± 0.15, and —3.03 ± 0.16, assuming the 
source distance of 300, 600, and 1,000 m, respectively. 



DISCUSSIONS 



With respect to the duration and its precedence to the lightning, the present event is 

n 

similar to one event observed by Torii et al. [151 ]. on another nuclear plant located on the 
same coastal area with a separation of 300 km. However, compared with their results, the 
present work provides a few new findings. That is, we have confirmed that the signal is 
dominated by 7-rays arriving from the sky direction. In addition, our results reveal that the 
photon spectrum extends up to 10 MeV. 

The photon spectrum detected in the present work, extending up to 10 MeV, clearly 
indicates a non-thermal emission mechanism. The most promising candidate process is 
Bremsstrahlung. Then, we can draw a following possible scenario to explain the present 
observation [171 ] . Generally, a winter thundercloud along the Sea of Japan is considered 
to have 3-layer (positive-negative-positive) charge structure in the vertical direction 19]. 
Then, the lower part of such thunderclouds is positively charged. Electrons, first produced 
by cosmic rays, will be accelerated by the runaway electron avalanche process [4J from the 
middle part to the cloud bottom (or top), emitting Bremsstrahlung 7-rays toward the ground 
(or sky). Since the 7-rays must be relativistically beamed into a narrow cone (~ 3° half-cone 
angle for a 10 MeV electron), we can detect such photons only under rare conditions when 
the cone happens to fall on our detector. 

Let us briefly evaluate energetics of the above scenario. Since we observed a fluence of 
1 — 10 MeV photons as n p = 2.4 x 10 4 m -2 , the total number of these photons at the source 
N p becomes <i 2 n p e _1 Af2, where d is the source distance, e is the atmospheric attenuation 
factor, and AQ is a solid angle of the possibly beamed emission. The values of e for 3 MeV 
photons are calculated, e.g., as 0.3, 0.08, and 0.01 for d = 300, 600, and 1000 m, respectively. 
Employing AQ = 8.6 x 10~ 3 (or equivalently a half-cone angle of ~ 3°) which is appropriate 
for 10 MeV electrons, we obtain iV p - 10 8 (rf/300m) 2 (0.3/e). Since a 10 MeV electron 
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traveling ~ 50 m in the atmosphere is expected to lose a minor fraction, 77 ~ 10~ 2 [20j, of its 
energy into 3 MeV bremsstrahlung photons, the overall energy of the produced relativistic 
electrons is estimated as 

E e = N p /r] x 10 MeV - 10" 2 (d/300m) 2 (0.3/e) J. (1) 

Of course, the estimate is subject to rather large uncertainties, because we may not have 
observed the peak of the radiation beam. Nevertheless, the estimated E e is a negligible 
fraction of an energy associated with a single lightning, 10 7 — 10 10 J {2]]], ensuring that the 
scenario is energetically feasible. 

In summary, we have successfully shown that electrons are accelerated beyond 10 MeV 
in thunderclouds, for an extended period, prior to lightning discharges rather than in co- 
incidence with them. Thus, this kind of observation of non-thermal x-/7-ray emissions 
associated with thunderclouds and lightning will provide a valuable key to the particle ac- 
celeration processes in strong electric fields. 
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of radiation safety group of Kashiwazaki-Kariwa power station, Tokyo Electric Power Com- 
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experiment. We are grateful to S. Otsuka and Y. Ikegami for design and construction of 
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FIG. 1: Schematic views of the two main sensors and the plastic scintillator of Detector- A. PMT 
means " photomultiplier" . Numbers represent scales in millimeter. 
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FIG. 2: Count rates per 20 seconds of 4 inorganic scintillators between 21:10 and 22:10 UT on 
2007 January 6. All abscissa represent universal time, and all errors are statistical la. Panels (a) 
and (b) show the count rates in > 40 keV data from the BGO and Nal scintillators of Detector-A, 
respectively. Panels (c) and (d) give > 40 keV and > 600 keV data from the Nal and Csl counters 
of Detector-B, respectively. The gaps in panels (a) and (b) are due to regular interruption of data 
acquisition of Detector-A every hour. 
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FIG. 3: Detailed count rate histories of the radiation sensors, and the light and electric field 
monitors, between 21:40 and 21:50 UT. All abscissa are universal time. All bins except those in 
panels (c) and (d) have la statistical errors, (a) Count histories per 6 seconds in > 3 MeV energies 
from the main Nal sensors of Detector-A, without (black) and with (red) the BGO anticoincidence, 
(b) 3 — 10 MeV counts every 6 seconds from the Nal (black) and Csl (red) scintillators of Detector- 
B. (c) One-second optical data variations, (d) One-second electric field data variations, (e) > 1 
MeV counts every 6 seconds from the plastic scintillator of Detector-A. 
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FIG. 4: Background-subtracted photon energy spectra accumulated over the 36 sec. All horizontal 
and vertical axes show the photon energy and the counts per each energy bin, respectively. Error 
bars are statistical la. (a) Spectra from Detector-A. Black and red circles show data of the main 
Nal counters (summed over the two units), without and with the BGO anticoincidence, (b) Spectra 
from Detector-B. Black and red squares are from Nal and Csl counters, respectively. Black and 
red dashed lines show predictions of the best-fit incident power-law model, determined over the 
energy range where the dashed lines are drawn. 
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